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We study the gravitational radiation reaction in compact binary systems composed of neutron 
stars with spin and huge magnetic dipole moments (magnetars). The magnetic dipole moments 
undergo a precessional motion about the respective spins. At sufficiently high values of the magnetic 
dipole moments, their interaction generates second post-Newtonian order contributions both to the 
equations of motion and to the gravitational radiation escaping the system. We parametrize the 
radial motion and average over a radial period in order to find the secular contributions to the 
energy and magnitude of the orbital angular momentum losses, in the generic case of eccentnc orbits. 
Similarly as for the spin-orbit, spin-spin, quadrupole-monopole interactions, here too we deduce the 
secular evolution of the relative orientations of the orbital angular momentum and spins. These 
equations, supplemented by the evolution equations for the angles characterizing the orientation of 
the dipole moments form a first order differential system, which is closed. The circular orbit limit 
of the energy loss agrees with loka and Taniguchi's earlier result. 



I. INTRODUCTION 



Neutron star and black hole binary systems are among the most probable sources of the gravitational radiation 
emitted in the frequency range of the Earth-based interferometric detectors such as the Laser Interferometric Grav- 
itational Wave Observatory (LIGO) H, VIRGO GEO 1^, and TAMA ^ (all will detect in the High-Frequency 
Band: 1 Hz to 10"* Hz) and the envisaged Laser Interferometer Space Antenna (LISA) 5, 6] . (Low- Frequency Band: 
lO^'* Hz to 1 Hz) respectively. A post-Newtonian (PN) treatment of 3.5 PN orders is generally agreed to describe 
with sufficient accuracy both the motion and the gravitational radiation up to the point of the innermost circular 
orbit (ICO), at least in the case of neutron star binaries jj. While the 3PN order approximation is expected to locate 
the ICO with an accuracy of 1% for binary systems with comparable masses, it breaks down in the region inner to 
the ICO in the latest stage of the final coalescence, when numerical relativity is required. However, there is a gap 
between the failure of the PN expansion and the berinning of the merger, this being called the intermediate binary 
black hole (IBBH) problem ^ , recently assessed in |3| . 

The effort of providing templates for the search of the gravitational waves in the noisy background is underway, at 
least for binaries on circular orbits, characterized by their mass and the radius of the orbit. The assumption of circular 
orbits is justified by the circularizing property of the gravitational waves 'lO'|, still interesting situations can occur in 
galactic nuclei, when freshly formed binary systems have not enough time to circularize ^11, 12]. The error produced 
by the neglection of the eccentricity during the search for gravitational waves was estimated to be substantial ^3 ■ 
The assumption of disregarding other characteristics of the binary components, like the spins, quadrupole moments 
and even magnetic dipole moments (to be discussed in the present paper), in many cases is not well founded. For 
example there is increasing observational evidence for an overhelming majority of galactic black holes (even colliding 
ones) with maximally allowed spins p^ . While neutron star binaries are not expected to be significantly spinning 
[TU l , 01 , neutron stars with super-strong magnetic fields (magnetars) have also been detected 0| . 

Therefore a generic treatment, allowing eccentric orbits and all types of contributions up to the 3.5*''PN order 
would be desirable in the PN-description of gravitational waves. The complexity of such results forbids their immedi- 
ate implementation in templates, however with increasing computational capacity this is hoped to become possible. 
An other complicating feature arising in the eccentric case is the occurence of several harmonics in the incoming grav- 
itational wave which, contrarely to the circular case (where the dominant frequency is twice of the orbital frequency) 
will have comparable intensities 0|. In the search for precessing binaries the use of effective templates has been 
proposed recently ^ , . 

The eccentricity was taken into account in the second PN order description of 21j. The 1.5 PN spin-orbit and 2PN 
spin-spin contributions to the orbital evolution and gravitational radiation reaction were considered in |2if23, 24, 2I 
I26ll27ll28ll29| . The gravitational quadrupole-monopole ty pe radiative contribution was also treated recently ,3^, the 
results matching the previous circular orbit computation |31j . 

The effect of possible strong magnetic fields of the compact binary components was discussed by loka and Taniguchi 
[3^ . for the case of circular orbits. They argue that in the upper limit of 10^^ G for the magnetic fields the coupling 
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of the magnetic dipole moments generates contributions of the same magnitude than the 2PN corrections to both the 
equation of motion and the gravitational radiation. They also compute the corresponding electromagnetic radiation, 
and find out that the electromagnetically radiated power is " much less" than the power radiated away gravitationally. 
The ratio of the two types of radiations is 3/64 cf. Eq. (16) of 's^l, which counts approximately half of a PN order. 
Being interested in the leading order radiation due to the magnetic dipoles, we will disregard the electromagnetic 
component of radiation. Whether the leading order gravitational radiation due to magnetic moments truely appears 
at second or only at higher orders depends on the strength of the magnetic fields. Actual observational evidence 
supports magnetic fields of 10^^ G for isolated magnetars ,17J. 

In the present paper we would like to extend our previous computations carried on for the spin-orbit (SO), spin-spin 
(SS) and quadrupole-monopole (QM) interactions for the case of the magnetic dipole - magnetic dipole interaction. 
This would mean a generalization of loka and Taniguchi's description allowing for eccentric orbits. 

In Sect. II we set up the formalism, by introducing the generalized true and eccentric anomaly parametrizations 
for the magnetic dipole - magnetic dipole perturbation. As proved earlier (33i | , all relevant integral expressions can be 
easily computed by use of the residue theorem, with the additional bonus that in the majority of cases the only pole 
is at the origin. Our treatment follows closely Refs. [2^, ^29] and 30J. In contrast with the energy the magnitude 
of the orbital angular momentum is not a conserved quantity in the absence of radiation at this order. Thus we 
introduce its angular average L in order to characterize the perturbed radial motion. The radial period, as well as 
the relation between L and the time average (L) resemble the results of the previously discussed cases. 

Sect. Ill contains the main results of the paper. These are the secular evolutions of E and L due to the magnetic 
dipole - magnetic dipole contribution to the gravitational radiation Also, as in previous cases, the evolution of the 
angle variables and 7, which characterize the relative orientation of the spin vectors and orbital angular momentum 
vector are derived. However the evolution of this set of variables [E, L, Ki, 7) does not close to a first order differential 
system, as in previous cases. In order to close the system, we have to compute the evolution equations for the angle 
variables Pi characterizing the relative orientation of the spins and magnetic dipole moments. The reliability of our 
results is checked in the Concluding Remarks, where in the circular limit we recover the expression for the energy loss 
given in [s^ . 

The velocity of light c and the gravitational constant G are kept in all expressions. 

II. THE RADIAL MOTION 

We consider a binary system composed of neutron stars with magnetic dipole moments dj (properly scaled in order 
to absorb all dimension-carrying constants). The dipole-dipole interaction is given by the Lagrangian |32j: 

C = Cn + Cdd = + + 1 [3(n • di)(n • da) - di • da] , (1) 

where v = f and r = rn. Accordingly, the total acceleration is the sum of the Newtonian a^v and the dipole-dipole 
contribution sldd: 

a. = ajy + sldd , 
Gm 

S.N = — , 

aoD = ^\ [di • d2 - 5(n • di)(n • d^)] n + V(n ■ di)dj I . (2) 

The orbital angular momentum L = r x p = Lat evolves due to the dipole-dipole perturbation as 

L= jE("-di)(nxdj) , (3) 

and its magnitude L changes accordingly 



i = LL=i^^(n.dO(v-fn).dj . (4) 



Lr- 



The energy E — rdC/dr — C and the total angular momentum J = L + Si + S2 are the constants of this motion. 
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We introduce three orthonormal coordinate systems /C and ]Q with the axes (c, L x c, L) and (bi, Si x bi, Si), where 
c and bi arc the unit vectors in the J x L and Si x L directions, respectively. In the system /C the polar angles and 
ipi of the spins arc defined as Si =(sin cos , sin Ki sin ipi , cos k^) (Fig 1) . We also introduce the polar angles and [3i 
of the the magnetic dipole moments di in the corresponding system /C* as di =(sinQ;i cos sin Ui sin [3i^ cos a^) (Fig 2). 
The transformation IQ ^ K represents a sequence of rotations Rz{—Ti)Rx{—Ki), where the angles = cos~^(c • bi) 
satisfy the relations + V'i = 7r/2 (see Fig. 2 for a proof). Thus in the K, system the vectors appearing in Eq. Q 
become^ : 



n 





di = di 



where we have introduced the shorthand notations 



Pi sin-0j + ai cos'ipi 
(Ji sin V'j - Pi cos V'i 




(5) 



Pi = sin cos /9i , 

Vi = sinaisin/3i , 

Ui = cos ai sin Ki + Ui cos Ki 

Q — cos OLi cos Ki — Vi sin Ki 

Since \s, independent of the velocity, = —C^dd holds: 

d\d2 



Edd — 



2r3 



[Ao-3B2ix)] 



(6) 



(7) 



^ Due to the magnetic dipole - magnetic dipole perturbation all three of the above defined coordinate systems may evolve (to DD-order). 
We do not consider these evolutions, as the components of the above vectors enter solely in various DD-terms, thus they are needed 
with Newtonian accuracy. 
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FIG. 2: The magnetic dipole moment of each neutron star undergoes a precessional motion about the spin vector of the 
respective neutron star. The projection of the magnetic dipole moment into the plane perpendicular to Si is denoted ^di, a 
vector characterized by the azimuthal angle /3i. On the intersection of the planes perpendicular to Ln and Si lies the vector 
bi. As the projection of Si into the plane of motion lies in the plane defined by Si and Ln, a plane to which the vector bi is 
perpendicular by definition, the relation Ti + tpi — n/2 holds. 



Here we have introduced the Newtonian true anomaly parameter x = ip — TpQ and 

^0 = 2cos A + 3(^1(72 - /02O-i) sin((5i - (52) - 3(/9i/92 + 0-10-2) cos((5i - ^2) , (8) 
^fe(x) = {cri(^2 - PiP2)cos{kx + Si + 62) - {pi<J2 + P2(^i)sm{kx + 61 + S2) , (9) 

where A is the angle subtended by the magnetic moments and Si = ipo — ipi. We stress here that he angular average 
of (O vanishes {Bk{x)) = for any integer k. 

With these notations, the evolution of the magnitude of angular orbital momentum is: 

L='-^B'2ix), (10) 

where a prime denotes derivative with respect to x- Iii DD expressions pi and ai can be regarded as constants, 
since the evolution of the angles Pi and Ki arc of 1.5PN order and ai are constants. Similar considerations hold for 
A. After the integration J Lx~^dx, we obtain 

L{X) - io-^^^{(3GmM + 4A)So 

-(3GmM + 4Acosx)B2(x) + AsinxB'^ix)} , (11) 
where A is the magnitude of the Laplace-Runge-Lenz vector for a Keplerian motion characterized by E and L 



A = ( G'm'p' + =^ ] . (12) 



From we see that L(0) = L{2tt) = Lq holds. Let L denote the angular average of L{x) 



L^Lo-^^^^iiGmfx + AA)Bo . (13) 
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Then we can express L{x) in terms of L rather than Lq 



L{x) ^ L + SLdd , (14) 
SLdd - ^^-:^[{3Gm^i + 4Acosx)B2ix)-AsmxB',{x)] • (15) 



2L3 

The magnitude of the velocity and the radial equation take the form 



v' - -[E-EnD{r,x)] + '^-^, (16) 
fi r 

= -[E-EDD{r,x)] + (17) 

Turning points are defined as solutions of the f = equation. They can be expressed as 

min 

where E^j^ = Ejjd (r±), 5U^j^ = SLdjj (r±) and r± are the turning points for a Keplerian orbit characterized by E 
and A. Substituting |7J and (|15() we obtain 

- + ^{iAT Gm,)A, + AB,] . (19) 

Next we introduce the generalized true and eccentric anomaly parametrizations of the radial motion, following the 
generic receipe from Ref. [83l |: 



/i(Gm^ + A cos x) 2AL2 (Gm^ + A cos x)2 
A = A[(3G2mV'+^')A + (G2mV'+^')Bo] 

+Gm^i [{G^m^fi^ + 3A'^)Ao + 2A^Bq] cos x , (21) 

r = ^^-^^ i + ^^{A(A + fio) + GmMoCose} . (22) 

Introducing suitable complex variables, radial expressions can be averaged by use of the residue theorem, as described 
in |33| . Following this recipe, the radial period turns out to have the Keplerian form 

\ 3/2 



T = 27rG™ (^-^ j . (23) 

(This holds whenever the radial dependence of the perturbing terms in the radial equation p7|l is either l/r^ or 1 /r^ 
[33 ■) Time average of L{x) gives 

{L)=L + \Gm^i\2G^m^^J? - SA^) - 2 {^^iiEf'^ L^j . (24) 

2/iyl^L3 L J 

This latter expression allows for an equivalent expression of all forthcoming results in terms of (L) rather than L, 
case needed. We remark that the square bracket in Eq. (|24l) coincides with the factor F1/F2 of the corresponding 
expressions in the quadrupole-monopole and spin-spin cases, Refs. (3fl| and p^ . 



III. LEADING ORDER MAGNETIC DIPOLE-MAGNETIC DIPOLE CONTRIBUTION TO THE 
EVOLUTIONS OF THE DYNAMICAL VARIABLES UNDER RADIATION REACTION 

A. Energy loss 

The radiative change in the energy to leading order is given by the quadrupole formula 
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where I^^^^^ is the 3'''' time derivative of the system's symmetric trace-free (STF) mass quadrupole moment tensor. 
To leading order it is given as 



^ fi{x^x'f^^ . (26) 

Inserting v'^ and r'^ from p6|) and (|17|l . together with Euu and L{x) from and 114|l into the quadrupole formula 
we obtain 

dE f dE\ f dE" 



dt \ dt J \ d.^ / DD 



(27) 



with the coefficients au given by 

01 = 3/iAr(-22G'TO/xV + 171^) , 

02 = 6(-llG2m2/zV+6£;ZV7-^ + 5GTOA*^L2r-5lZ^) , 
aa = -^lAr{22Gm^?r + SIL^) , 

a4 = 2Z2(-6£;^r2 -5Gm/iV + 39L^) . (30) 

To compute the averaged loss of the energy we parametrize the above expression by the true anomaly x and pass to 
the complex variable z — exp(ix). The only pole of the integrand is located in the origin. Using the residue theorem 
we obtain 

'dE\ _ /dE\ /dE\ 



dt I \ dt / \ dt I ' 
dE \ G-'vfi{-2E^fl^ [USE^L^ + 732GV + 425GV V) , (32) 



dt / 15c5L7 

where the coefficients Gi^2 are 

G2 = 708E^L^ + 1QQ20G^ E^L^m^fi^ + 18865G'^EL^m'^ fi^ + 8316G^mV^ • (34) 



B. Change in the magnitude of orbital momentum 

Since there is no secular spin evolution in the 2PN order, the loss in the magnitude L under radiation reaction can 
be written as 

dL - dJ , , 

— ~ L • — , 35) 

dt dt ' ^ ' 

where ~ denotes equality modulo spin terms, which average out due to {dSi/dt) — 0, see j29j |. The instantaneous loss 
of the total angular momentum J to leading order is 

ClL oc 

Inserting w^, and the respective components of n and di from jSJ into H35|) we obtain 

dt \ dt) \ dt) ^jj' 

3 1 



dt 
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with the coefBcients given by 

bi = 3fiAr{2GEmn^r^ - EL^^ir^ - 9GL^m^i^r + 8L^) , 

&3 = ^Ar{2GEmii^r^ + SEL^^it^ - 9Gl?m[i^r - 2U^) , 

bi ^ L^{~l8Efj.r^ -8Gmf/r + 15L^) . (38) 
After averaging, we obtain the secular loss in the magnitude of orbital angular momentum: 

'dL\ /dL\ /dL\ ^^g^ 

dt/, ' 

dL\ Gdid.2(-2Enf/^ r . n 

where the coefheients D12 are given as 

Di = -6fiA'^{31EP + 90G'^m'^fj.^) , 

L>2 = 2b2E'^L^ + l2QQG'^EUm'^^? + 8Q^G'^m^^i^ . (42) 



/ \ dt / \ / jjj^ 



dL\ 4G2m(-2i;^)3/2 



C. Evolution of the angles characterizing the spins under radiation reaction 

The relative orientation of the momenta can be described by the angles Ki and 7, see Fig 1. Consequence of 
(dSi/dt) = 0, the angle 7 is conserved: 

^cos7~0. (43) 

The angles were found to evolve due to both the spin-orbit and the spin-spin interactions . Here we compute 
a third contribution to their evolution, due to magnetic dipole-magnetic dipole interaction: 

— COSKj ~ ^(Si-LcOSKj)- — 

dt J DD L \dt J OD 

^ ^ sin Ki ^ { X! "fc ['^J cos(A:x + 5^ + 5^) - sin(fcx + 5, + 5j)] 



^ j=i k=i 



-I-M4 sinx [cTj sin((5i — <5j) — pj cos{5i — 5j)] + U5 [pj sin((5i — Sj) + (jj cos((5i — <^j)]| (44) 



with the coefficients Uk 



ui = -U3 = pAr{2Epr^ - 3i^) , 
U2 = -2L^{Gmfi^r + 3L^) , 
U4 = 2p.Ar{2Epr^ - SZ^) , 

M5 = 2L^{4:Epr^ + Gmp^r - SZ^) . (45) 
Averaging yields the following secular expression for the change of the angle k^: 

/dn,\ 3Gd,d2{-2Ef,)y^ ^ WA^AM 



\dt I 5c5Z9 



-I- V2 [cTj cos{6i -Sj) + Pj sin((5j - Sj)] } , (46) 

where Vi^2 are 

Vi = bpA^iAEL^ + 7G^m^^r') , 

V2 = 48£;^Z* -t- 140G^£;Z2toV^ + TOG^m^fj.^ . (47) 

Note that although the detailed expressions of the secular evolutions of E, L and Kj are different from the corresponding 
expressions charcterizing the gravitational quadrupole-monopole interaction j33|, the coefficients ai_4, 6i_4, Ui_5, 
C'1,2, £'1,2 and Vi,2 are identical! This is related to the similar structure of the respective Lagrangians. 
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D. Evolution of the angles characterizing the magnetic moments under radiation reaction 

The first order differential equations (|^. and togetlier with the algebraic constraints (presented in 

detail in [l^]): 

Si sin Ki cos V'l + <S'2 sin K2 cos ■(/'2 = , (48) 
cos Ki cos K2 + sin Ki sin K2 cos (?/'2 ^ V'l) — C0S7 . (49) 

does not form a closed system, due to the presence of the angles A, a^, and (3i contained in Aq, Bq, and Q. This is 
a new feature to be contrasted with the spin-orbit, spin-spin and quadrupole-monopole cases. Therefore we need the 
radiative evolution equations for the above angles either. 

For this purpose we remark that A can be expressed in term of the other enlisted angles as 

cos A = {p20-i ~ pi(T2) sin((5i - S2) + {P1P2 + (^1(^2) cos((5i - 82) + C1C2 ■ (50) 

A second remark is that similarly to (dSi/dt) — 0, we expect (ddi/dt) = to hold. The assumption can be lifted by 
considering any neutron star model which relates the magnetic dipole moment to other characteristics, like the spin^ . 
Therefore the radiative change of the angles = cos~^(Si • dj) is beyond the order of accuracy of our computation 



dt 



cos tti ~ 



(51) 



and the only task remains to compute the radiative evolution of 

-Ldi • fSi X L 



cos Pi 



sm Ki 



where by 



sin (3i cos Ki cos -iAi + cos (3i sin tpi 
"'"dj = I sin Pi cos Ki sin ipi — cos Pi cos V'i 
— sin Pi sin Ki 



(52) 



(53) 



we have denoted the unit vector aligned with the projection of the magnetic dipole moment to the plane perpendicular 
to the spin of the respective neutron star. The vector product of the projection of the magnetic dipole moment and 
the direction of the spin, appearing in the subsequent expressions, is readily obtained 



sin Pi sin -0^ — cos Pi cos t/i^ cos 
"di X Si = ( — sin Pi cos V'i — cos Pi sin V'i cos Ki 

cos Pi sin Ki 



(54) 



Up to terms which do not contribute to the secular expressions 

d 



dt 



-""di X Si dJ cos Pi dL cos Ki „ dKi 

cosPi ~ , , . — - -rr^-jT - -■ cos/3» — 

L [X) sm Ki dt L [X) dt sm Ki dt 



By virtue of Eq. H35|l and Eq. (|52(l we find that to the leading order in the DD-terms, Eq. (|55|l becomes: 



d 
di 



{sin Ki cos Pi) ~ {^^' ^ ^' 



-di X S 



llL}( 



dJ 

K'dt 



(55) 



(56) 



DD 



As {d3/dt)N ^ L (for which the expression in brackets acts as a projector) we find that similarly to the angles Ki, the 
radiative evolution of the angles Pi receives no Newtonian contribution. Employing Eq. I|44(l . the detailed expression 



^ Actual magnetar models relate the non-radiative spin precession to the magnetic dipole-moment. Co ncerning the evolution of the 
magnetic field, they deal onl y wi th non-radiative evolution, ranging from allowing no evolution at all l34l to a non-linear magnetic field 
decay through the Hall-drift l35l . 
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of Ea. H56|l turns out to be 



2 



d ^ 3Gdld2 SinPi-^ (yr-^ r • /, r cm 

Tt ^ - 5^5^2X2^7 1^ 1 2^ + + '^J ) + ^"^('^^^ + '5. + 5, )] 

— U4 sin X [fj cos((5i — + Pj sin(6i — 6j)] + [aj sin((5i — Sj) — pj cos(^i — 6j)] | (57) 



and its average over a radial period gives 



/dPi\ 3Gdirf2(-2S/i)-'^/2 2 

(1),,- Sc^Z^sin^. [p,cos(^. + ^,)+a,sm(^. + ^,)] 

+ V2 [cj sin(5i — Sj) — Pj cos((5i — Sj)]] . (58) 

Note that the coefficients appearing in the instantaneous and averaged loss of the angles Pi agree with the coefficients 
H45|l and H47|) in the radiative change of k^. 

IV. CONCLUDING REMARKS 

We have given a closed system of differential equations governing the evolution of a set of dynamical and geometrical 
variables under the magnetic dipole - magnetic dipole contribution to the radiation reaction. These are the energy, 
magnitude of orbital angular momentum, angles between the spins and orbital angular momentum and angles between 
the dipole moments and the respective spins. These evolutions add to previously derived PN, 2PN, SO, SS and QM 
contributions. For sufficiently strong magnetic fields this new contribution is of second PN order. 

As a check of the above results we compare the circular orbit limit of the energy loss with the previous result of loka 
and Taniguchi 32]. Imposing the circularity condition in an average sense (as described in Ref. '30|), the following 
relations hold 

En =E^Edd = -^P^ , L% = = Gm^i^ro . (59) 

Here tq is the radius of the unperturbed circular orbit and the angular average Edu of O is 

did2Ao 

3 
'^0 



Edd = . (60) 



Inserting E and L from Eqs. H59|l into the expression of the energy loss H31|) we obtain the radiative loss of energy 
for the previously defined circular orbit: 

dE\_ 32GWm' a 12rfirf2.4o \ 

dt / ScVg Gm/irg ) ' ' 

Comparison with the corresponding result of 32] is achieved by remarking that Eq. (jS)) can be written in the form 

A = 3(L-di)(L.d2)-di.d2 , (62) 

and by computing the averaged radius of the quasi-circular orbit in terms of the radius of the Newtonian circular 
orbit tq: 

After the required series expansion we find that our energy loss, when specialized to circular orbits and expressed 
in terms of (r), is in perfect agreement with the energy loss contribution due to magnetic dipole - magnetic dipole 
interaction, given in Eq.(f6) of Ref. (3^. 
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